Breakdown of the initial value formulation of scalar-tensor gravity and its physical 

meaning 
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We revisit singularities of two distinct kinds in the Cauchy problem of general scalar-tensor 
theories of gravity (previously discussed in the literature), and of metric and Palatini f(R) gravity, 
in both their Jordan and Einstein frame representations. Examples and toy models are used to shed 
light onto the problem and it is shown that, contrary to common lore, the two conformal frames are 
equivalent with respect to the initial value problem. 
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I. INTRODUCTION 

The 1998 discovery of the acceleration of the cosmic 
expansion, obtained by studying type la supernovae [H, 
spurred an enormous amount of activity on dark energy 
models, mostly based on cosmological scalar fields. Cer- 
tain models are set in the context of scalar-tensor gravity 
instead of Einstein's theory, and are dubbed "extended 
quintessence" 0. Moreover, as an alternative to postu- 
lating a mysterious form of dark energy, various authors 
( [!, E| > see @ f° r a review) have considered the possibil- 
ity that the acceleration of the universe is caused instead 
by a modification of gravity at large scales: the Einstein- 
Hilbcrt action 



Seh = 



2 k 



d A x 



is generalized to 
1 



S 2k 



d i x^f(R) + S^[g ab ,^} 



(1) 



(2) 



where f(R) is an arbitrary, twice differentiable, function 
of R. Here k = 8irG, G is Newton's constant (that 
will be unity, together with the speed of light, in the 
geometrized units employed), R is the Ricci curvature, 
S (m) = j d A x^g~C( m) [g ab ,4>] is the matter part of the 
action, ip collectively denotes the matter fields, and we 
follow the notations of Q . 



If the action @ is varied with respect to the metric g ab , 
one obtains the metric formalism with fourth order field 
equations [1, 5 ; if the metric and the connection T% c are 
considered as independent variables (i.e., the connection 
is not the metric connection of g a b), but the matter part 
of the action does not depend explicitly on T, one 
obtains the Palatini formalism with second order field 
equations If, instead, depends on T, one obtains 
metric-affine gravity 

It has been shown [§] that metric f(R) gravity is dy- 
namically equivalent to a Brans-Dicke (BD) theory with 
BD parameter ujq = 0, while Palatini f(R) gravity is 
equivalent to an ujq = —3/2 BD theory. The general form 
of the scalar-tensor action, of which BD theory 1(J 11 1 
is the prototype, is [l2T ] 



Sst = J d 

+ SW[g ab ,i>] 



f(<t>)R "(0 



V c 0V c 0- V(0) 



(3) 



where is the BD-like scalar field and /(</>) > is re- 
quired in order for the effective gravitational coupling 
to be positive and the graviton to carry positive kinetic 
energy and not being a ghost. V((f>) is the scalar field po- 
tential, while f((f>) and to((f>) are two (a priori arbitrary) 
coupling functions. BD theory is recovered as the special 
case /(</>) = (f> and Lo((f>) = wo/4>, with loq =const. The 
field equations derived from the action §5§ are 



/(<£) ( Rab - ~ 9ab R) = cj(<f>) ( V a 0V fc - ^ 5ab V c 0V c ) - Vg ab + V Q V b / - g ab Uf + T { ™ ] , 



3(f 



2/ 



□</>- 



J 3/'/" , uf \ f 2Vf 

^ + —+v r r c4>V(t> =y T+2V - — 



(4) 
(5) 



r 



where a prime denotes differentiation with respect to 0, □ = .g ab V a Vfc, and T ab = —j2= 



2 



The original motivation for BD theory was the im- 
plementation in relativistic gravity of the Mach princi- 
ple, which is not fully embodied in general relativity, by 
promoting Newton's constant to the role of a dynami- 
cal field determined by the environment fiol . IeO ] . Later 
on, it was discovered that string theories and supergrav- 
ity contain BD-like scalars: in fact, the low energy limit 
of the bosonic string theory (which, although unphysi- 
cal because it does not contain fermions and is not su- 
persymmetric, was one of the early string theories) is 
indeed an wo = -1 BD theory [13}. Moreover, BD the- 
ory can be derived from higher-dimensional Kaluza-Klein 
theory, higher dimensionality being an essential feature 
of all modern high energy theories. A p-brane model in D 
dimensions leads, after compactification, to a BD theory 
with parameter [14| 



tion: the field equation for the BD scalar is 



w 



(D-l)(p-l)-(p+l) 2 
(U_2)(p-l)-(p + l)s 



(6) 



These properties have renewed the interest in BD and 
scalar-tensor gravity since the 1980s, following the rise of 
string theory. However, the more recent surge of interest 
in scalar-tensor gravity that we are witnessing is moti- 
vated by cosmology and is linked to attempts to explain 
the present cosmic acceleration (see 15j, ll6| for reviews 
of scalar-tensor gravity in the cosmology of the early and 
present universe). 

Motivated by the past and recent interest and also 
by developments in numerical relativity, the initial value 
problem of scalar-tensor gravity was studied by Sal- 
gado [17[ who, using a first order hyperbolicity analysis, 
showed that the Cauchy problem is well-formulated for 
theories of the form 



J> - V(cf>) 



S (m) [cjabM 



(7) 



when 5( m ) is "reasonable" 51 1 and is well-posed in vacuo. 
It was then straightforward to generalize this work to 
BD theories with constant BD parameter luq ^ 1 which, 
in turn, was used to show that the Cauchy problem of 
metric f(R) gravity (equivalent to an loq = BD the- 
ory) is well-formulated and well-posed in vacuo, while 
the Cauchy problem for Palatini f(R) gravity (equiva- 
lent to an u>o — —3/2 BD theory) is not well- formulated, 
nor well-posed [lj|. A second paper by Salgado and co- 
workers using a second order hyperbolicity analysis [20] 
showed the well-posedness of to = 1 theories, and the 
extension to to =const. theories (with the exception of 
to = —3/2) is straightforward because the principal part 
of the field equations does not depend on to. 

In retrospect, it is easy to see why the u>o = —3/2 BD 
theory does not admit a well-posed initial value formula- 



te) 



V 



2V 



(8) 



and reduces to a first order constraint when luq — > —3/2. 
Technically, this fact prevents the substitution of 
back into the equations for the other dynamical variables 
in order to eliminate second derivatives of 4> and spoils 
the reduction to a first-order system (see [l9[ for details). 
In practice, the second order dynamical equation for the 
variable 4> is lost when uio = —3/2, <f> then plays the role 
of a non-dynamical auxiliary field and can be assigned 
arbitrarily a priori. Uniqueness of the solutions is then 
lost, as infinitely many prescriptions for <j> correspond to 
the same set of initial data. 

The present paper serves various purposes. First 
(Sec. II), we revisit the Cauchy problem for BD the- 
ories and, in particular, for the equivalent of Palatini 
f(R) gravity by using a completely independent approach 
based on the transformation to the Einstein conformal 
frame. This approach was deliberately avoided in pre- 
' J. While interesting in itself as 



17, m 



vious papers 

an independent check of previous results, this approach 
has the additional merit of fully establishing the physical 
equivalence between Jordan and Einstein frames at the 
classical level. These conformal frames have been shown 
to be equivalent in various other respects, and it would 
only make sense that their equivalence extend to the 
Cauchy problem. However, there are explicit statements 
in the literature, and much unwritten folklore, pointing 
to the contrary. We show here that the two frames are 
indeed equivalent, which removes previous doubts and 
fully establishes equivalence at the classical level; how- 
ever, this does not guarantee physical equivalence at the 
quantum level 2l|, [22 1 . 

The main purpose of this paper, however, consists of 
the study of the Cauchy problem for scalar-tensor the- 
ories of the general form Q and of two distinct types 
of singularities that may appear in their field equations. 
These theories were not covered explicitly in previous lit- 
erature, although the extension of the results of 17, 2(3] 
to include them is relatively straightforward. In addition, 
it is handy to consider the general form © of the theory 
in order to specialize the results to any scalar-tensor the- 
ory simply by prescribing specific forms of the coupling 
functions /(</>) and to{(f>) and of the potential V(<f>). We 
approach the problem in both the Jordan frame (Sec. Ill) 
and the Einstein frame (Sec. IV) obtaining, of course, the 
same results. 

In general scalar-tensor theories, there are two kinds of 
singularities to deal with: those at which f(4>) = 0, and 

a second kind identified by fi(4>) = to(4>) + 3 ^f/(ff = ^' 
which generalizes the loq = —3/2 pathology encountered 
in BD theory and in Palatini f(R) gravity. Singularities 
of the first kind should normally be excluded by requiring 



3 



that f(4>) > for all values of (j>, and this requirement is 
sometimes made explicit in the general formalism (e.g., 
[23j]); nevertheless, works incorporating these singulari- 
ties recur often in the literature, especially in cosmology. 

At the singularities of the second kind fi — (which 
have been known for a long time in particular incarna- 
tions of scalar-tensor gravity 0, [2^, [2(| ) , a phenomenol- 
ogy similar to that of Palatini f(R) gravity spoils the 
Cauchy problem for special forms of the coupling func- 
tion oj(4>), or for critical field values. While the scalar 
field is allowed to pass through these "singularities" 
in an isotropic Friedmann-Lemaitre-Robertson- Walker 
(FLRW) universe [U 0, [26[ , the points where = 

are known to give rise to curvature and shear singularities 
in the anisotropic case. These singularities were discov- 
ered in the special case of nonminimally coupled scalar 
field cosmology (corresponding to /(<ft) = j , — £,4> 2 and 
u> = 1) in the early universe 24, 2^, 27, 28^ and also in 
black hole perturbations [2^]. They also appear in the 
search for exact wormhole solutions with nonminimally 
coupled scalar fields [3C| ■ Singularities of both kinds were 
discussed in a more general context in [31], [32|. After 
clarifying and further generalizing this situation from the 
point of view of the initial value formulation in Sec. Ill, 
in Sec. IV we revisit this subject in the Einstein frame, 
exposing a situation analogous to uq = —3/2 BD the- 
ory (this is not merely an analogy, since the latter is a 
special case of the former). Finally, in Sec. V, we study 
nonminimally coupled scalar field theory as an example, 
recovering certain known properties and placing them in 
a general context. Sec. VI and VII contain illustrative 
toy models and the conclusions, respectively. 



II. EINSTEIN FRAME DESCRIPTION OF 
BRANS-DICKE AND PALATINI f(R) GRAVITY 

In this section we recall the definition of Einstein con- 
formal frame and show explicitly the non-dynamical role 
of the scalar field in the Einstein frame representation of 
the scalar-tensor version of Palatini f(R) gravity. This 
is necessary as a first step to understand the more in- 
volved situation that we will be facing in later sections 
with general scalar-tensor theories of the form ([3]) . 

The conformal transformation 



Jab 



9ab 



g a b 



(9) 



into its Einstein frame representation 

R 1 



Sbd = J d x y<-g 
£ (m) [n^gab,^]' 



2k 2 



where 



u{4>) 



<f> 2 (4>) 



u{4>) 

(12) 
(13) 



and a tilde denotes rescaled (Einstein frame) quanti- 
ties. The scalar field redefinition (fTQ|) breaks down when 
u)q = —3/2 and the scalar field <fi then remains undefined. 
However, eq. ^ still holds and one can write the Ein- 
stein frame version of the BD equivalent of Palatini f(R) 
gravity as [f| 



~> Palatini 



2 k 



d 4 x 



R 



v{4>) 



S (m) [<p~ 



(14) 



and the scalar field redefinition <f> — ► <f> with 



using the variables (g a b,4>). In this action, the scalar 
field <p does not play any dynamical role: it only acts 
as a factor rescaling the metric in but it has no 

dynamics, does not couple to R, and no kinetic energy 
of <f> appears in (|14jl . <p can be assigned arbitrarily in 
infinitely many ways not governed by the usual second 
order differential equation and, therefore, uniqueness of 
the solutions is lost. On the contrary, for any value of the 
BD parameter u> ^ —3/2 (in particular for the ujq = 
equivalent of metric f(R) gravity), the action is reduced 
to (TT2"]) , which describes a scalar field 4> coupling mini- 
mally to the curvature and nonminimally to matter. In 
vacuo (£( m ) = 0), this coupling to matter disappears and 
we are left with the action of Einstein gravity plus a min- 
imally coupled scalar field with canonical kinetic energy: 
it is well-known that this system has a well-posed initial 
value formulation 0, [l8| . The non- vacuum case is con- 
sidered later in Sec. V as a special case of more general 
scalar-tensor theories. 



III. THE CAUCHY PROBLEM FOR GENERAL 
SCALAR-TENSOR THEORIES IN THE JORDAN 
FRAME 



|2w - 



2 k 



map the Jordan frame action of BD theory 
Sbd = / d 4 x yf^g 



(10) 



(11) 



Let us now restrict ourselves to the Jordan frame and 
consider general scalar-tensor theories described by the 
action ([3]). These can be reduced to the action 



S 



d A x 



— i- 1 V C $V C $ - [/(*) 



<PR 
2 



_ s (rn) 
(15) 
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containing a single coupling function u*(<f>) by setting 
$ ee JO), w*($) = and [/($) = The 

actions ([3]) and (|15[) are equivalent if f(4>) is invertible 
with regular inverse but this does not happen if 

f'{4>) vanishes somewhere. 

It is also possible to recast BD theory as one in which 
the kinetic term of the scalar field is canonical, i.e., with 
uj = 1. Beginning with the action (fT5|) and setting 

$ = F((p) , (16) 

where the function F(ip) is defined by the equation 

F(<p) 



2(f 



(17) 



(fT5]) can be rewritten as 



S 



d x y/—g 



F(ip)R 1 



+ S (m) [ffo&,^] 



(18) 



where W(</?) = y [F(i^)]. This alternative form of the BD 
action can not be obtained when F(ip) does not admit a 
regular inverse F~ x (e.g., when dF/dip = 0). This is the 
case, for example, when F(tp) is represented by a series 
of even powers of (p [H, HI] . Note that (fT5)) is the form of 
the action considered in the studies of the Cauchy prob- 
lem 17, 2(J. In what follows, to achieve full generality, 
we discuss the action with two coupling functions. 
Moreover, there are two types of "singularities" to con- 
sider: we introduce them here and we will refer to them 
for the rest of this paper. In addition, one must distin- 
guish betwen two very different situations: that in which 
these "singularities" occur in an entire four-dimensional 
domain of spacetime, and that in which they occur only 
on hypersurfaces. Morever, we will approach all of the 
above from the two viewpoints of Jordan frame and Ein- 
stein frame. 



Singularities of the first kind 

The first type of singularities is identified by f(4>*) = 
and occurs for critical values 0* of the scalar field (if so- 
lutions of this equation exist). This equation may be 
satisfied in an entire four-dimensional spacetime region, 
or on a hypersurface. At a first glance, the former case 
seems rather trivial: in fact, naively, the effective gravita- 
tional coupling read off the action (f2"6"]l is G e // = l//($0- 
However, a more careful analysis of the effective gravita- 
tional coupling in a Cavendish experiment, which is the 
only one directly accessible to local experiments, yields 

13 



GeffW 



2ujf + {2df/d(j>f 



8tt/ 2w/ + 3 



(19) 



This expression can also be obtained from the study of 
cosmological perturbations [36] ] . The first type of singu- 
larities /(</>*) = corresponds to diverging effective cou- 
pling G e ff(4>) and separates regions in which G e // has 
opposite signs describing attractive or repulsive gravity, 
respectively. Stated this way, it may seem nonsensical 
to consider such values 4>* of the scalar. For example, 
it looks plain silly to consider, in BD teory, a spacetime 
region in which = 0, which makes the term c/)R/2 dis- 
appear from the BD action and corresponds to infinite 
strength of gravity. Nevertheless, there are examples in 
which exact (and non-unique) solutions of the field equa- 
tions have been found with (f> constant and precisely equal 
to in a region, or in the entire spacetime manifold 



(see [26J, |37J, |38( for examples in cosmology and [30| for 
wormhole solutions). Are these to be discarded a pri- 
ori! Perhaps not, because what is clearly unphysical are 
regions in which G e // < and the graviton is a ghost. 
Although rather pathological, regions in which G e f / is di- 
vergent may still be interesting in exotic situations when 
the birth of the universe or the interior of a wormhole 
are considered. Furthermore, these regions may still be 
relevant from the mathematical point of view if one is 
interested in finding exact solutions that, as simplified 
toy models, exhibit particular properties of scalar-tensor 
gravity. 

Let us come now to the more interesting situation in 
which f(4>) = on an hypersurface. This situation seems 
more reasonable, however such hypersurfaces separate re- 
gions of attractive from regions of repulsive gravity; in 
the latter, the graviton carries negative kinetic energy, 
a physically unacceptable property [23|, |39| . This fact 
seems to be forgotten in scalar-tensor theories more gen- 
eral than BD theory and with more freedom in the form 
of the functions f(4>) and to((f>). Papers in which G e // 
is negative or infinite have appeared surprisingly often 



ppcarea surpi 

I m, H m 



41, m 



over the past thirty years [2£ 
sometimes, such critical hypersurfaces <f>* are approached 
asymptotically [H^]. 

Let us proceed, for the moment, by adopting a purely 
mathematical point of view in the consideration of the 
Cauchy problem. When f(<f)) = 0, eq. ((3]) for the metric 
tensor degenerates. At these spacetime points the trace 
of eq. (|4|) becomes 



3\3<j) +(uj + if") V c 0V c + W - T = 



(20) 



Substitution of the value of obtained from this equa- 
tion into the second field equation ((5|) yields 



2 

7 



V c 0V c 



+ ^p(T-W)-2V'} 



(21) 



Knowledge of the values of <fi and of its gradient V c (j> 
on the hypersurface / = determines the Ricci curva- 



5 



ture. However, the equation for R ao disappears there, 
which means that all metrics with the same value of R 
satisfy the (degenerate) field equations on this hypersur- 
face: uniqueness of the solutions is lost and this surface is 
a Cauchy horizon. The initial value problem breaks down 
at these critical points. Therefore, even if we decide to 
allow the unphysical region G e // < by attempting to 
propagate initial data given in a G e // > region, we en- 
counter a hypersurface on which G e // — > oo which acts 
as a barrier and the initial value formulation ceases to be 
well-posed. 



Singularities of the second kind 

Let us introduce now a second type of critical values 
of the scalar field that have previously been associated 
to physical (curvature) singularities and that also corre- 
spond to a breakdown of the initial value problem. Fol- 
lowing the lesson of u> — —3/2 theory [19j], one notices 
that disappears from the field equation ([5]) when 



h{4>) = ^{4>) 



2/(0) 



= 



(22) 



Again, one has to distinguish two cases: a) eq. (|2"2"|) is sat- 
isfied in a four-dimensional spacetime region, and b) it 
is satisfied on a hypersurface. The former corresponds 
to regarding eq. ([22]) as specifying a particular form of 
the coupling function uj((f>) (given f{4>)), while the latter 
corresponds to seeing eq. (|2"2"j) as a trascendental (or alge- 
braic, depending on the forms of the functions lj and /) 
equation that may admit as roots special critical values 
<j) c of the scalar (f> [53| • 

Let us consider case a) first: this is completely analo- 
gous to the case of u) = —3/2 BD theory which eq. (j2"2"|) 
generalizes. When fi(<fi) vanishes identically for all val- 
ues of </> in a four-dimensional spacetime domain, the 
dynamics of the scalar 4> are lost together with □</> and 
with the second order of the partial differential equation 
for (/). The exception consists of situations in which the 
scalar satisfies □</> = 0, in which case there may be non- 
trivial dynamics for </>, but this quantity disappears spon- 
taneously from the field equations for the other variables. 
This situation includes general relativity with <\> =const. 
(for which the initial value problem is well-posed [6fl and 
the previous discussion obviously does not apply), and 
harmonic 0- waves. 

Situation a) is, of course, the only possibility when to 
represents a constant parameter instead of a function, as 
in BD theory. The general scalar-tensor theory is richer 
and allows one to contemplate the possibility b) that 
eq. (|22|) is satisfied on a hypersurface. It is interesting 
that, in the absence of matter, invariants of the Riemann 
tensor diverge at this hypersurface for anisotropic met- 
rics, while no such divergence occurs in isotropic FLRW 



spaces [25j, l3l| . Mathematically speaking, if /((/>) =/= 
and fi((/)) is a continuous function, a hypersurface where 
fi{4>) — separates two regions corresponding to oppo- 
site signs of fi (unless the form of fi is pathologically 
fine-tuned): in each of these, the Cauchy problem may 
be well-posed but when one tries to propagate initial data 
through such a hypersurface, D0 given by 

f /lj' 3/'/" , 



f'T , 2V , 2Vf 



2/ 



/ 



3(f) 2 
2/ 



(23) 



diverges. We have, therefore, a Cauchy horizon that is 
not hidden inside an apparent horizon, as in black holes, 
and where the theory crashes. The two regions separated 
by the hypersurface fi{4>) = are, again, disconnected 
by a singularity in the gravitational coupling G e f / (0) . 

To summarize this section: when the coupling func- 
tions /(</>) and u}{4>) are such that f((f>) = or f\{<j>) = 0, 
the initial value formulation breaks down and either the 
theory is unphysical because <f> becomes a non-dynamical 
auxiliary field, or the hypersurface /i(</>) = is a Cauchy 
horizon. In the first case, the problems found for Palatini 
f(R) gravity in [45| re-surface. The situation in which 
eq. (|2"2")l is identically satisfied is the generalization to ar- 
bitrary scalar-tensor theories of the situation already seen 
in oj = —3/2 BD theory and in Palatini f(R) gravity. The 
trace equation ([20f allows one to replace the trace T with 
an expression containing second derivatives of 4>. Then, 
the metric depends on derivatives of the scalar field of 
order higher than second and discontinuities, or irregu- 
larities, are not smoothed out by an integral of matter 
fields giving the metric g ao (for example, as in the usual 
Green function integral in the weak-field limit), but they 
cause step-function discontinuities in the metric deriva- 
tives and curvature singularities where the same matter 
distribution in Einstein's theory would generate a per- 
fectly regular geometry. 



IV. GENERAL SCALAR-TENSOR THEORIES 
AND THE CAUCHY PROBLEM IN THE 
EINSTEIN FRAME 

We now examine the initial value problem of general 
scalar-tensor gravity in the Einstein frame. The confor- 
mal transformation 



g ab -> g ab = ft 2 g ab , ft - y/fW (24) 

and the scalar field redefinition 



(25) 
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bring the Jordan frame action 



5 = / d 4 x \/—g 
into its Einstein frame representation 



(26) 



S 



d 4 x 



1 



C (r. 



f 2 



R 

2k 2 ■ 

(27) 

where U(</>) = V{^))/f 2 and / = f &(</>)). Again, 
apart from the now familiar coupling of the "new" scalar 
<f> to matter described by £( m ) // 2 (with the exception of 
conformally invariant matter), this action describes gen- 
eral relativity with a canonical scalar field which couples 
minimally to the curvature but nonminimally to mat- 
ter. As before, it is clear that the system has a well- 
posed initial value formulation in vacuo. This conclusion 

applies where the Einstein frame variables (g a b, (fij are 

well-defined, i.e., for /(</>) ^ and fx(<f>) 0. It can 
be shown that the Cauchy problem is well-posed in the 
presence of matter as well: this was already pointed out 
in ref. 23], but is checked at the end of this section by 
extending the first order hyperbolicity analysis of [l7j |. 

The exception is when fx((f>) — 0, in which case the 
scalar <f> can not be defined using eq. (|25|) . In this case, 

one can use the variables (g a b, </>) instead of (jj a b, 4>) , ob- 
taining the Einstein frame action 



R_ 

2k 



V(<t>) 



/ 2 W>) f 2 (<P) 



(28) 



Again, there are no dynamics for <f> and the Cauchy prob- 
lem is not well-formulated, nor well-posed, in this case 
due to the loss of uniqueness of the solutions. Morever, 
this result holds in both the Jordan and the Einstein 
frames, which then become physically equivalent in this 
respect. 

The breakdown of the scalar field redefinition (f25f is 
accompanied by other signals that something is going 
wrong with the physics when f\{4>) = 0. The effective 
gravitational coupling (p~9j) diverges when /i = (as a 
special case, G e // = ^o+3W> diverges as w o — * —3/2 in 
BD theory). Moreover, it changes sign when <j> crosses a 
critical value <j>* or <fi c . These critical values are attained 
by the scalar field in certain early universe inflationary 
scenarios with nonminimally coupled scalar fields, cor- 
responding to /((/>) = - — £<j> 2 and u> = 1 (£ being 
a dimensionless coupling constant) when < £ < 1/6 
24, 25, 26|, 2q |. The same phenomenon in more general 
scalar-tensor theories is considered in 

The authors of 



31, 3 



31[ find that, in Bianchi cosmolo- 



gies, the regions of the phase space at which fi(4>) — 
correspond to geometric singularities with divergent 
Kretschmann scalar R a bcdR ahcd ■ The /i = singular- 
ity is dynamically forbidden in a closed or critically open 



FLRW universe under the assumptions p > 0, V(<f>) > 0, 
and u > 1 [13|. 



The lesson of [25J, |28|, |3l| is that, if there is even a 



small anisotropy, the change from attractive to repulsive 
gravity at fx — can only occur through a shear or cur- 
vature singularity which stops the evolution of the geom- 
etry: nature's message seems to be that gravity can not 
spontaneously become repulsive in the absence of exotic 
matter violating the energy conditions (it is the purely 
gravitational sector of the theory that we are studyng 
here). 

Note that in the theories considered by [U [12] , which 
have u> = 1, the singularity f\ = is automatically re- 
moved by requiring that f(4>) > (i.e., that the graviton 
is not a ghost); however, this is no longer true when the- 
ories with ui not identically equal to unity are considered 
and critical values <p c of the second kind can still occur 
even when f((f>) > V0 — but this necessarily requires 
u < 0. 

Let us now extend the first order hyperbolicity analysis 
of [ljj to Einstein frame scalar-tensor gravity. We follow 
closely, and adopt the notations of, [lj, [H)} in order to 
facilitate comparison, setting k = 1. The Einstein frame 
field equations are 



Gab = V Q ^V b 



1 

-j9ab" 



: ^V C 0- U(4>)g a b 



T 



(m) 



i 2 m)) 



fab{4>] +f 



(m) 
ab 



dU(4>) 
d<f> 







(29) 



(30) 



Because T a b[0] does not contain second derivatives of <j>, it 
is possible to give a first order formulation as in general 
relativity. The nonminimal coupling factor 1 / f 2 ! (<^(<j6)) 
multiplying on the right hand side of eq. ([2^]) does 
not generate derivatives of 4> and therefore is immaterial. 

The 3 + 1 ADM formulation of the theory defines the 
usual lap se, shift, extrinsic curvature, and gradients of 
4> d, [l?! . Assuming the existence of a time function t 
such that the spacetime (M,g a b) is foliated by a family of 
hypersurfaces S t of constant t with unit timelike normal 
fi a , the 3-metric is defined by h a b — g a b + n a fib and /i" 
is the projection operator on St. The relations n a h a = 
— 1, h a bfi b = h a bn a — 0, and h a hb c = h ac are satisfied. 
Further introducing the lapse N, shift vector iV a , and 
spatial metric hij , the metric is written as 

ds 2 = - [N 2 - N l N^ dt 2 - 2N i dtdx i + h./iyd.r' (31) 



(i,j = 1,2,3), with N > 0, n a = -NV a t and 

N a = -hlt b , (32) 
where the time flow vector t a satisfies t a V a t — 1 and 

t a = -N a + Nn a (33) 
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so that N — —h a t a and N a n a = 0. The extrinsic curva- 
ture of Yi t is 

Kab = -h c a h d b V c fi d . (34) 
The 3D covariant derivative of h ao on S 4 is defined as 



A w T ai - fcl ... 



h ai 



(35) 



for any 3-tensor (3) T ai - fcl ... , with Dji ab = 0. The spa- 
tial gradient of the scalar field and its momentum are 



and 



respectively, and 



Q a = -D a (j 



n = Cn<j) = h c \7 c (j) , 



Ki 



1 / dhi 



2N \ dt 



i>..\, • n,s, , 



ti=l(d t 
N V 



d t Q t + N l diQ z + Qid z N l = A (iVfi) . 
The stress-energy tensor is 3 + 1-decomposed as 

Tab = Sab + Jaflb + Jb^a + Eh a hb , 

where 

S a b = h a hfT c d — S ab [4>] + Slffl , 

J a = ~h C a f cd h d = J a {0}+J { a m) , 

E = n a n b f ab = E[4>] + , 



(36) 
(37) 

(38) 

(39) 
(40) 

(41) 

(42) 
(43) 
(44) 



and f = S — E, where T is the trace of T a b and S is the 
trace of S ao . The Gauss-Codacci equations provide the 
Einstein equations projected tangentially and orthogo- 
nally to T, t as the Hamiltonian constraint @, 17 1 



(45) 



{3) R + K 2 l\,J\ ,J = 2E , 
the vector (or momentum) constraint 

Dik\ - D.K = Ji , (46) 
and the dynamical equations 

d t k) + N l diK) + KidjN 1 ~ K^diN 1 + D'DjN 



N 



- {3) mN - kkk) - 



S - E ) 5) - 2S) 



,(47) 



where k = k\. The trace of this equation yields 

dtK+N^ik+^AN-NKijk^ = — (§ + e) , (48) 

where = DWi. 

Further introducing Q = Q C Q C , one computes 

E[$\ = ~ (n 2 + Q 2 ) + U(4>) , (49) 

J[4>] = -IlQa , (50) 



\ (q 2 n 2 ) + u&) 



, (51) 



SafcM = QaQb ~ h a b 

while 

S[$\ = \ (sri 2 - Q 2 ) - 3U(4>) (52) 



and 



S[$\ - E[4>] = (n 2 - Q 2 ) - MJ{$) . (53) 



The "total" quantities entering the right hand side of the 
3 + 1 field equations are then 



1 ~ 9 l~ 9 



E = -Q' z + -U z + U(^) + E^ 



Ja = -flQa + , 



(54) 



(55) 



S a b — h a b 



\ (q 2 - n 2 ) + 1/(0) 



while 



QaQb + Sif , 



s = -m$) f - ^- + s {% 



S-E = 

fl 2 - Q 2 - AU{4>) + S {m) - E {m) , 

S + E = 2Y1 2 - 2U(<f>) + S {m) + E {m} . 

The Hamiltonian constraint becomes 

fr 2 n 2 
2 2 

= E^ + U{4>) , 



(56) 



(57) 



(58) 



(59) 



(60) 
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while the momentum constraint (|4G|) is 

bik\ - bik + flQi = J| m) , (6i) 

the dynamical equation (|47|) is written as 



3t_E + N l diK) + kid 3 N l - kfaW + frDjN 



-^RjN - NkK] + ^-2U{^)5) + JVQ'Q 



_ N 
~ ~2 

with trace 



/ 2S 3 



(m) z 



(62) 



9 t AT + N l diK + (3) A TV - NKtfK^ 



-ATT 



N 



(-217(0) + S' (m) + £ (m) ) (63) 



where 17 



£ fi f[ - ft£ - Q C I) C (lniv) - Z) C Q C = -□<£ 



(64) 



In vacuo, the initial data (hij , Kij ,<j>,Qi, flj on an initial 
hypersurface Eg obey the constraints ((60|) and (|61jl plus 



Qi - Di<l> = , 



(65) 



In the presence of matter, the variables 



m) 



T(m) 
<J a , 



and 



£ab are also assigned on the initial hypersurface. Fixing 
a gauge corresponds to prescribing lapse and shift. The 
system contains only first-order derivatives in 

both space and time once the d'Alembertian D<f> is writ- 
ten in terms of </>, V c 0V c </>, and their derivatives by using 
cq. ([6411 . From this point on, everything proceeds as in 
Ref. [171 ] and the nonminimal coupling factor /(</>(</>)) in 
f^f =Ti; n) // 2 does not have consequences because it 
contains no derivatives of si™ , 



ah 



or E( m >. The re- 
duction to a first-order system indicates that the Cauchy 
problem is well-posed in vacuo and well-formulated in 
the presence of those forms of matter for which it is well- 
formulated in general relativity. We do not duplicate Sal- 
gado's analysis here, and we refer the reader to 17, 13] 
for details. 



Equivalence between conformal frames 

At this point, it is clear that the initial value formula- 
tion is well-posed in the Einstein frame if it is well-posed 
in the Jordan frame, and vice-versa. The two frames are 
equivalent also from the point of view of the Cauchy prob- 
lem, contrary to folklore and recurring statements in the 



literature. To this regard, it is often remarked that the 
mixing of the spin two and spin zero degrees of freedom 
g a b and <fi in the Jordan frame makes these variables an 
inconvenient set for formulating the initial value problem, 
which is consequently not well-posed in the Jordan frame, 

while the Einstein frame variables (jj a b , </>^ admit a well- 
posed Cauchy problem completely similar to that of gen- 
eral relativity. (A rather casual remark in the well-known 
paper [3^] (see also the more recent Ref. [23} ) seems to 
have been quite influential in this respect, without fur- 
ther questioning of it in later literature until the recent 
work of Salgado 17J). In the light of this work, which 



is carried out completely in the Jordan frame, the stan- 
dard lore is obviously false. Old works also hinted to the 
fact that the Cauchy problem is well-posed in the Jordan 
frame for two special scalar-tensor theories: Brans-Dicke 
theory with a free scalar <p [3] , and the theory of a scalar 
field conformally coupled to the Ricci curvature [47j . The 
implementation, in the Jordan frame, of a full 3+1 for- 
mulation a la York [48j for use in numerical applications 
further dispels the myth that the Cauchy problem is not 
well-posed in the Jordan frame • 

Were this folklore true, the Jordan and Einstein frames 
would be physically inequivalent with regard to the 
Cauchy problem, but we have shown that this is not the 
case. In fact, the equivalence between the two conformal 
frames does not break down even when the scalar field re- 
definition <j> — > 4> fails. The Jordan and Einstein frame are 
still equivalent, with respect to the initial value formu- 
lation, for general scalar-tensor theories and, therefore, 
they are equivalent at the classical level, thus dissipat- 
ing residual doubts left in this regard in 22j. However, 
the two conformal frames seem to be inequivalent at the 
quantum level ([5l], [i^] and references therein). 



EXAMPLE: THE NON-MINIMALLY 
COUPLED SCALAR FIELD 



We arc finally ready to consider, as an example, the 
theory of a scalar field coupled nonminimally to the Ricci 
curvature. In fact this example, many features of which 
are well-known, has sometimes already guided us through 
this paper. The action is 



Snmc = I d 4 x y/^g 



-V{J>) 



Mm) 



(66) 



where £ is a dimensionless coupling constant (in our no- 
tations, conformal coupling corresponds to £ = 1/6), and 
a m is a suitable coupling constant. The field equations 
are 
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(1 - <0 2 ) Gab = K 

Oj> - %■ - £Rcj> = 



V a 0V & - -g a bV c (j)V c (j} - V{4>)g ab + £ (&,&□ - V a V 6 ) (0 2 ) + 



(m) 



(67) 
(68) 



(see 15|, |49[ for a discussion of alternative ways of writing 
the field equations). By neglecting the matter part of the 
action, the Ricci curvature can be eliminated from the 
Klein-Gordon equation obtaining 



1 + (6j - 1) <0 2 
1 - <0 2 



1 - n^ 2 



:i-6ov c </»v c 



(69) 



Singularities of the first kind correspond to f(cf>) = i 
£0 2 = 0, or to the critical scalar field values 



±1 



(70) 



and can only occur if £ > 0. They correspond to diverg- 
ing effective gravitational coupling 



G. 



G 



eff 



1 - K £(j) 2 ' 



(71) 



which changes sign if the scalar crosses ±0*. The re- 
quirement f((f>) > avoids these critical values. How- 
ever, one could decide to momentarily ignore the physi- 
cal interpretation of the theory and to allow these critical 
values from a purely mathematical point of view; then, 
the latter return to haunt the Cauchy problem and pre- 
dictability. 

The quantity f± is, in this theory, 



h{4>) = 



l + < (6g- l)cj) 2 

1 - K^ 2 



(72) 



The roots of the equation fx = 0, which exist if < £ < 
1/6, are the critical values of the second kind 



±1 



^(1-60 



(73) 



The non-uniqueness of the solutions and the breakdown 
of the Cauchy problem marked by the critical values 
±0»,±</> c are seen as follows. When <fi = <Po =const. 
and matter is absent, the theory reduces to vacuum gen- 
eral relativity with a cosmological constant, and the field 
equations reduce to 



G a b + kg ab = 



A 



i + <0 2 



(74) 



(75) 



The trace of eq. ([74]) gives R ■■ 
eq. (|75p in turn implies that 

R 



4 A which, compared with 



£00 



(76) 



If <fr = ±0, (the critical values of the first kind), then 
it must be Vb = and, therefore R a t = 0. The 
Klein-Gordon equation yields the extra necessary con- 
dition Vq = 0. All vacuum solutions of general relativity 
(Rab — 0) are also solutions of the field equations ([6"7| 
and (HHJ) with = ±0* . 

If instead = ±0 C (the second kind of critical values), 
eqs. J7JJ) and ([75]) yield 

K (1 - 60 V(±<t> c ) 



and 



A = 



V\±4> a ) = T 



2(1 - 30 



2V<(1-6QV(±^ C 
l-3£ 



(77) 



(78) 



At these critical scalar field values of the second kind, 
the dynamical equation ([69]) for </> loses all the second 
derivatives of <j> (contained in Dip) and, consequently, the 
dynamics for this field (except for special solutions satis- 
fying □(/> = 0). In isotropic FLRW spaces, solutions are 
known which cross the critical values ±0 C , or cf> is identi- 
cally equal to one of these values. However, the situation 
can be worse: there are physical curvature and shear sin- 
gularities in anisotropic Bianchi models 25, 28, 31, 50l ]. 
Moreover, Barcelo and Visser [3(| find diverging Ricci 
scalar R for spherically symmetric wormholc solutions. 
These examples correspond to solutions which cannot 
cross the barrier (j) = ±<p c . 

The conformal transformation to the Einstein frame is 
9ab — > gab = ^ 2 g a b with 



n = yj\ - ni4> 2 



(79) 



and the redefinition bringing the scalar field into canon- 
ical form is 



dej) 



y/l-Kg(l-6Q 
1 - <0 2 



(80) 



By integrating the last equation, the Einstein frame 
scalar cj> can be explicitly expressed in terms of <j) as 



£V6«0 2 + yfi - £ (1 - 6g«g 

£V6«0 2 - v 71 - £ (1 - 6£) k</> 2 



f{4>) , 

(81) 
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where 



Let us first consider the point particle action 



,, , 1 - en 1/2 

f{4>)= [ K £ i 



(^(1-60 n^) (82) 



for0<(< 1/6 and 



/WO 



6£- 1 



for £ > 1/6, while 



1/2 



arcsinh 



(vm^iy^) (83) 



for £=1/6. 

The Einstein frame action is 



if M> 



(85) 



n[<t>)^ m) } (86) 



where 



U 







r 



0(0) 



i - <0 2 (0) 



(87) 



and 



1 - <0 2 (0) 



2 • 



When <p = ±0*, the conformal transformation of the met- 
ric breaks down, while the redefinition of the scalar field 
becomes invalid when <f> — ±0 C . In this last situation, one 
can still use the variables (g a b, <ft) to define an Einstein 
frame in which the action is simply 



S = I d A x 




v(4>) 



(1 - K£pf (1 - K^f 



C (m) 

(89)' 

with no dynamics for <p, which becomes an auxiliary field 
and can be assigned arbitrarily. 



VI. TOY MODELS 

In this section we consider toy models in the context 
of point particle dynamics, which help obtaining some 
insight into the "singularities" of the first and second 
kind of scalar-tensor theories. 



5' 



dtL(x(t),x(t),y(t),y(t)) 



dt 



i 2 f(y) w(y)y 2 



-J(x) 



(90) 



where an overdot denote differentiation with respect to 
the time t, the generalized coordinates x and y mimic 
the metric g a b and the scalar respectively, the func- 
tions f(y) and w(y) represent /(</>) and u>{4>), while J 
represents the matter sources. Since we are interested in 
the purely gravitational sector, we will set J to zero in 
most of the following. 

The coordinate x is cyclic and the Euler-Lagrange 



equations | (f^ 



|£=0(i = l,2) yield 



±f{y) = c 



(91) 



, y , w'{y)y 2 , f'{y)x 2 , . 

w{y)y + — g — + — 2 — = ' ^ ' 

where a prime now denotes differentiation with respect 
to y and C is an arbitrary integration constant. 
i) The analogue of a singularity of the first kind f(4>) — 
in a domain is f(y) = on an interval, which implies 
C = and 



, „ w'(y)i/ 2 

Hy)y + = o 



This equation admits the first integral 
rv 



dy' VWl/Tl = Ci(t- t ) 



(93) 



(94) 



where C\ and to are integration constants. Let us con- 
sider now, for the sake of illustration, the choice w(y) = y 
yielding the solution 



y(t) = C 2 (t-U) 2/3 , 



(95) 



with Ci and i* integration constants. Note that, because 
/ = 0, there is no equation for x(t) and the dynamics of 
this variable are lost: the initial value problem is not 
well-posed because x(t) can be assigned arbitrarily and 
is not determined uniquely by initial data (xq,Xq) at an 
initial time to. 

ii) Let us consider now the situation in which f(y) van- 
ishes at isolated points mimicking the critical scalar 
field values 0*. Then, in the system (f9"Tj) and (|9"2"|) , ei- 
ther C = or else x — > oo as y — > y„. If C ^ 0, then 
x = C/f(y) — > ±oo as f(y) — + ± and, therefore, also 
x(t) — ► ±oo or the solution is not of class C 1 and its 
derivative docs not exist. In the first barrier sep- 

arates the regions f(y)>0 and f(y) < 0, however special 
solutions which traverse the barrier y = y* can in princi- 
ple exist. 
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If C — 0, an exceptional solution x(t) =const., y(t) = 
Cz(t — t*) 2 ^ 3 passes through this barrier, however this 
corresponds to the special value C = and it disappears 
when C ^ 0. 

Hi) Let us consider now the case w(y) = on an interval, 
corresponding to a singularity of the second kind fx (y) = 
on a domain. Then, we are left with 



*/'(!/) = 



(96) 



(97) 



From eq. (|97|). either x(t) =const. and then it must be 
C = with no equation left to determine y(t), or the 
equation f'(y) = is an algebraic (or trascendental, but 
not a differential) equation that determines constant val- 
ues y* of y (if it admits roots). Assuming that f(y*) ^ 0, 



then x*(t) 



f(v.) 



t + x . The solutions (x*(i), y*(t)), if 



they exist, are the only ones and correspond to excep- 
tional initial conditions and, in this sense, there are no 
dynamics for y. 

iv) We can now consider the situation in which w(y) van- 
ishes at isolated points y c , mimicking isolated singulari- 
ties of the second kind fi(<p c ) = 0. Consider, for exam- 
ple, the choice w(y) = y, f(y) = y — 1, for which the 
system (|9H and reduces to 



(98) 



x(y-l) = C, 



•9 -2 
V X 

yy + y + (y - 1) y = o 



(99) 



Assuming that y is not identically unity, it is yy 



2(y-i) 



= 0; at y = one has y c = ±C and one can not 
assign arbitrary initial conditions on the "hypersurface" 
analogue y = 0, but only the initial data (xq, Xq, yo, yo) = 
(x , — C, 0, ±C) are allowed there, where C and xo are 
arbitrary constants. The region allowed to the dynamics 
in the four-dimensional space (xo, io, yo, yo) is only two- 
dimensional, due to the presence of the first integral (|9"Tj) 
and of the additional first integral 54 1 



w(y)y 2 



C 2 



= H = const. 



(100) 



If a solution attains the critical value y = 0, it must as- 
sume the values (xo, xo, yo, yo) = i x o, —C, 0, ±C) there, 
for which the "energy" Ti. can only take the values 



H 



C 2 



[w(0) + 1] 



(101) 



(where eq. (HJll) has been used) everywhere along the or- 
bits of the solutions. 



The analogue of the Einstein frame 

Let us consider again the toy model action (|90|) ; the 
tranformation to the Einstein frame for scalar-tensor 
gravity is a change of variables modelled by the trans- 
formation (x,y) — > (£,77) defined by 



= Vf¥)dx 



dr) = y/w{y) dy 



(102) 



(103) 



In terms of these new variables, the action (|90[) is rewrit- 
ten in the canonical form 



5' 



dt 



e 
2 



(104) 



which mimicks the Einstein frame representation of the 
scalar-tensor action with the matter sources J now de- 
pending on both the "new metric" £ and the "scalar field" 
(y, or 7/ through y(r))). Zeros of either f(y) or w(y) make 
the analog of the conformal transformation plus scalar 
field redefinition (fT02|) . (TT03| ill-defined. Moreover, if 
f(y) 7^ an( i onr y w (y) vanishes, one can still consider 
an "Einstein frame" representation with the variables £ 
and y, in terms of which the action is simply 



S 



dt 



i 2 

j-J&y) 



(105) 



It is clear that, similar to the case considered before for 
scalar-tensor gravity, there are no dynamics for the vari- 
able y, which can be assigned arbitrarily [55| . This is bad 
news if this variable plays a physical role because there 
are no equations to rule it and it can only be assigned 
from outside the theory, which is akin to invoking a mira- 
cle to produce any effect that one may desire and results 
in a complete loss of predictive power for the theory 



Singular points of ODEs 

To conclude this section, we comment on the fact that, 
in the theory of ordinary differential equations (ODEs), 
it is rather common to encounter situations in which the 
phase space is divided into two disconnected regions, with 
only exceptional solutions, or a restricted submanifold 
of solutions, crossing the boundary between these two 
regions. Consider, for example, the ODE 



t 2 y 



2y = 



(106) 



which has t = as a regular singular point. Two linearly 
independent solutions are 



ft® 



t 2 



V2{t) 



(107) 
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The first solution crosses undisturbed the t = "barrier" , 
while the second cannot (that is, t — is a barrier to 
at least some of the solutions). Consider also the third 
solution in (— oo, 0) U (0, +oo) 



V3(*) = 



if t < , 



t 2 if t > . 



(108) 



2/3 is continuous with its first derivative at t = (but 
the second derivative is not defined there). Now, y%(t) 
and 2/3 (t) are linearly independent solutions which sat- 
isfy the same initial conditions (2/(0), 2/(0)) = (0,0) at 
t = 0. These two otherwise distinct solutions intersect 
at the origin of the phase space, which signals the break- 
down of the initial value formulation at t = 0. It is not 
surprising, therefore, that for the more complicated sys- 
tems of partial differential equations ruling scalar-tensor 
theories, the Cauchy problem breaks down at the ana- 
logue of singular points of the equations. Depending on 
the particular form of the coupling functions f{4>) and 
ui ((/)), special solutions crossing the barrier may or may 
not exist. The situation in ODE theory in which no such 
special solution exists is exemplified by the equation 



t y + 5ty + 3y = 



(109) 



the physics becomes unpredictable. Physically, this is as- 
sociated to a change in sign of the effective gravitational 
coupling (|19[) . which diverges at both kinds of critical 
points. It seems that nature abhors such changes from 
attractive to repulsive gravity (and vice-versa) which, 
formally, only take place through a singularity of G e f / . 
This, however, says nothing about exotic forms of matter 
which can source repulsive gravity through the field equa- 
tions, a completely different and seemingly perfectly le- 
gitimate mechanism from the mathematical point of view 
(although the violation of all the energy conditions would 
certainly be questionable on physical grounds). 

To conclude, we remark that a possible cure for the 
problem of Palatini f(R) gravity (already outlined in 
Refs. [HM [5(1) could be the insertion into the gravita- 
tional action of terms that introduce higher order deriva- 
tives into the field equations. Then, the dropping out of 
□</> from the field equations will be immaterial. However, 
unless such higher derivative terms appear in the Gauss- 
Bonnet combination, they will introduce ghost fields. 
A study of the initial value problem for these Gauss- 
Bonnet-corrected theories will be presented elsewhere. 

VF acknowledges Marcelo Salgado for a useful com- 
ment. This work is supported by the Natural Sciences 
and Engineering Research Council of Canada. 



For t > and for / < 0, two linearly independent solu- 
tions are yi(t) — l/t and 2/2 (t) = l/t 3 . No choice of the 
arbitrary constants C\ t 2 in the general solution 



y(t) 



C\ C2 



(110) 



in (—oo,0) U (0, +00) produces a solution crossing the 
barrier t = 0. 



VII. CONCLUSIONS 

In principle, two kinds of "singularities" for the Cauchy 
problem are possible in general scalar-tensor theories: 
those ("first kind") at which f(<p) = 0, and those 
("second kind") at which fi(<j>) — 0. Although state- 
ments that these should be rejected outright have been 
voiced in the literature [2j|, solutions corresponding to 
critical values of the BD-like scalar field of both first 
H M, S M, and second kind @, [U, HI, EL Hfll 
have been studied in the literature. Critical points of 
the second kind may appear benign when studied in 
a spatially homogeneous and isotropic FLRW universe, 
but they reveal their true nature of geometrical sin- 
gularities when analyzed in anisotropic Bianchi models 
[25I [2R l3ll [Hoj . Here, following recent developments in 
the theory of the Cauchy problem of scalar-tensor grav- 
ity, we have shown that the latter is not well-posed at any 
of those critical points. The solutions are not unique and 
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